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NOMENCLATURE 


a^ = Semi-major axis of the satellite before capture maneuver 

aj = Semi-major axis of the satellite after capture maneuver 

a^ = Semi-major axis of the 'dead’ body before capture maneuver 

a£ = Semi-major axis of the 'dead' body after capture maneuver 

c t = Constants ; i = 1,2, 6 

e^ = Eccentricity of the satellite before capture maneuver 

e| = Eccentricity of the satellite after capture maneuver 

e^ = Eccentricity of the 'dead' body before capture maneuver 

e£ = Eccentricity of the 'dead’ body after capture maneuver 

E = Eccentric anomaly 

G = Universal Gravitation Constant 

H = Constant ; Eqn < Z . 3 > 

A A 

i,j,k = Unit vectors along 77 , £ axes , respect 1 ve 1 y 

k = Kinetic energy factor 

k. = Constants ; 1 = 1,2 4 

K = earth's gravitational constant 

1^ = Length of the tether between c.m. and m^ 

1^ = Length of the tether between c.m. and m 2 

1^. = Length of the tether 

m = Mass of the hook 

o 

rrij = Mass of the satellite 

m 2 = Mass of the 'dead' body 

M = Mean anomaly 

n = Angular velocity 
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c.m. of the two massive particles (Fig. 4) 
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Constants ; l = 1.2 j Eqns. (2.1) - (2.2) 

c.m. of the tethered system (Fig. 4) 

the radius vector of c.m. of the tethered bodies 

the radial distance ; l = 1,2 (Fig. 4) 

the radial distance of the hook from the earth's 

centre (Fig. 1) 

the radius vector of the satellite from the 
earth's centre. (Fig. 2) 

the radius vector of the 'dead' body from the 
earth's centre. (Fig. 2) 

the radius vector ; i = 1,2 ; ) = 1,2 (Fig. 4) 

1 1 me 

time of perigee passage ; l = 1,2 
components of the vector u i ; l = 1,2 

velocity of the satellite at the time of capture maneuver 
velocity of the ’dead' body at the time of capture maneuver 
unit vector along r ; l = 1.2 

radial and transverse unit vectors, respectively 
unit vectors along x, y axes, respectively 
velocity of the satellite after capture maneuver 
velocity of the ’dead’ body after capture maneuver 
unit vector along 1 ^ ; l = 1,2 

components of the vector v ; l = 1,2 
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angular separation between and c.m. <Fig. 2) 

_ A 

oi^ = the angle between r. amd v^ <Fjg. 4> 

ft = angular separation between m^ and c.m. (Fig. 2> 

_ _ ^ /S t*. 

<5r = perturbations in r, <5r = 6% 1 + 6r) j + <5( k 

AAA 

= <5x u + <5y u + 6 z k 
x y 

e = constant 

<p - angle between r and v^ (Fig. 4) 

y = direction of v ; i = 1,2 

l l 

X = constant ; l = 1,2 in Eqn (2.2a> 

/j = mass of the massive particle P ^ 

v = mass of tethered bodies ; i = 1,2 

i 

v = constants in Eqn (3.7a), (3.4a), (3.4), respectively 

yj = pitching angle at time ' t ' (Fig. 2, Fig. 4) 

yj Q - pitching angle at time t = 0 

yj = constant 

© = angle between OP and r (Fig. 2) 

© ^ = angle between OP and r^ (Fig. 1) 

= angle between OP and r^ (Fig. 2) 

© 2 Q = polar angle of the 'dead' space object as measured from 

the tethered satellite position at t=0, the time of 
perigee passage 

©^ = angle between OP and hook (Fig. 1) 

dot = derivative with respect to time 

prime = derivative with respect to the true-anomaly © 
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ABSTRACT 


Two important dynamic aspects of tethered satellite 
system (TSS) are considered here. In the first part, the 
attention is focussed on the extraction of kinetic energy from the 
fast moving space wandering objects to the relatively slow moving 
satellites. The study clearly establishes the usefulness of small 
tethers that enable a significant momentum transfer from space 
debris to the satellite. 

Finally, the dynamics of the two-body TSS with its c.m. 
at one of the stable Lagrange points namely or L,. in the 
gravitational fields of the two massive bodies is investigated. 
The results obtained here show that in general the orbital 
perturbation due to finite tether length causes the c.m. of the 
tether to slowly drift away from these stable points. The 
perturbations are only of the second order, however. The orbital 
drifts over even relatively large life-spans of the TSS turn out 
to be insignificant. 
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CHAPTER - 1 


INTRODUCTION 


1.1 Preliminary Remarks : 

The advent of Space Shuttle extended the possibility of 
exploration and exploitation of the space activities. Tethered 
Satellite System <TSS) has been conceived as one of the devices to 
extend the capability of space shuttle to perform scientific/ 
applications investigations and opera t ional activities. The 
concept educes a multiuse tethered system with closed-loop 
control, capable of supporting a payload or satellite suspended 
from the shuttle cargobay toward or away from the earth and 
distances upto approximately 100km from the shuttle. These 
satellite payloads, although tethered to the shuttle, can be 
placed in their new orbits permanently, retrieved, refurbished and 
reconfigured for subsequent missions. A detailed description and 
operation of these has been presented by Rupp and Laue^\ 

The successful missions which involved TSS were carried 

( 2 > 

out aboard the Gemini spacecrafts . In these missions, the 
manned Gemini XI and Gemini XII space vehicles were tethered to 
the unmanned Agena Vehicle. The former used the rotating 
configuration that was maintained by Gemini thruster reactor 
control system and the latter used gravity gradient stabilized 
conf i gurat i on . 

In 1990 it was proposed to fly first tether 
depl oyer/retr i eval system in space. But this mission delayed due 
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to temporary setback to the shuttle programme that followed the 

<3 4) 

Challenger disaster. The first of _these missions ' , developed 
and funded by NASA and Italian National Space Plan (CNR) is now 
proposed for the year 1992-93. It will deploy approximately a 300 
Kg subsatellite upward from the shuttle on a 20mm thick kevlar 
tether of 20km length, where the gravity gradient will stabilize 
the two masses along a vertical configuration. It will 
investigate the interaction of tether - which is designed to have 
a conducting core - with the earth's magnetic and ionospheric 
environment. A second mission will be flown about two years 
later, using the same deployer system aboard the shuttle at 220km 
altitude deploying a subsatellite 100km downward to investigate 
earth’s atmosphere. Several other tether or tether related 
applications are in incipient stages. The 1990's should see 
atleast half a dozen such missions which will demonstrate specific 
tether capabilities and increase our understanding of their 
benefits and behaviour in space^\ 

Umpteen number of possible tether applications have been 
proposed for a variety of purposes. Some of the applications are 
briefly summarized below: 

1 . Electromagnetic Generator : 

An insulated conducting tether connected to a space 
craft and possibly terminated with a subsatellite can generate DC 
electrical power to supply primary power on board. The DC 
electrical power is generated by induction through motion of the 
conducting tether in the geomagnetic field. 
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2. Space station platform with controlled gravity laboratory : 


A linear, three-mass-conste 1 lat 1 on that consists of the 
space station and a science/technology platform at opposite ends 
of a tether in a gravity gradient stabilized configuration is 
proposed. A platform on an elevator or crawler mechanism would be 
made to move along the tether length between its two ends for 
servicing. Alternately, the platform can be used as controlled 
gra i 1 ty laboratory wherein the g-level can be regulated by moving 
the elevator. 

3. ULF /ELF /VLF Communications Antenna : 

The variations in current through insulated conducting 
tether connected to a spacecraft and terminated at both the ends 
with plasma contactors can produce ULF/ELF/VLF waves 
enabling worldwide communications. This tether antenna can be 
self-powered using the current induced in it by power or 
externally powered through an on-board transmitter. 

4. Satellite Boost from Orbiter : 

A satellite is deployed along the tether line upwards 
(i.e., away from the earth) from the shuttle orbiter. Momentum 
exchange takes place from shuttle orbiter to the satellite. Now 
if the tether is severed, the satellite will move to a higher 
orbit and at the same time the shuttle will be deboosted somewhat 
before it is finally brought back to the earth. 
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5. Tether Rendezvous System: 


The tether rendezvous system will be using a 'Smart' 
hook which can capture and retrieve payloads, space shuttle to 
space station. 

1.2 A Brief Review of Literature : 

A tether deployed from a space station is attached to a 
docking module. This module can capture and retrieve the shuttle 
allowing a remote rendezvous. This shuttle-docking by tether has 
been discussed by Hunter^^ . 

A decayed satellite can be either repaired or reboosted 
into a higher orbit by a rendezvous maneuver with the help of a 
permanent tether attached to a shuttle. This eliminates the need 
for launch of a new satellite. The above rendezvous maneuver is 
also discussed by Hunter^*. 

The idea of a tethered satellite approach for attaining 
low Earth orbit was put forth by Smithsonian Astrophys 1 ca 1 
Observatory in 1974^*. In their proposal to the Atmospheric, 
Magnetospher i c and Plasmas in Space <AMPS> science working group, 
they described the use of long tether system called 'Skyhook' for 
low altitude research. 

A tether rendezvous system can capture and retrieve 

payloads, OTV's or space shuttle to space station with the help of 

( 8 ) 

a 'Smart' hook by a rendezvous maneuver. This tether 
rendezvous system helps in saving of a great amount of propellant 
for incoming vehicles. 
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' By using long tethers, transportation of mass from 

higher or lower orbits can be achieved. The concept of skyhook 
tether transportation system has been discussed by Penzo^\ 

In many-body problem, knowing the initial positions and 

velocities of three or more massive particles moving under their 

mutual gravitational forces, one can determine their final 

positions and velocities at any time. This problem was first 

< 9 > 

formulated by Newton . The particular solutions to the 
Three-body problem were given by Lagrange in 1772***** . His 
circular restricted three-body problem deals with how the two 
massive particles moving in circular orbits about their common 
centre of mass affect the motion of the third infinitesimal 
particle under the influence of their gravitational fields. 
Poincare, Hill and others**** made a great deal of study of the 
circular restricted 3-body problem, 

1.3 Purpose and Scope of Investigation : 

The investigation proposed here deals with two important 
aspects related to a tethered satellite system. First, the 
possibility of using tether for the extraction of kinetic energy 
<KE) from the fast moving space-wandering 'dead’ objects to the 
relatively slower satellites involved in useful space applications 
has been studied. Given the initial conditions - namely the 
positions and velocities of the two bodies - it is proposed to 
find out the final velocities of the two bodies. This study 
assumes that the two bodies are moving in the same plane. The 
satellite is tethered to a small hook through which it captures 
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the 'dead' spacecraft passing by. The tethered system is assumed 
to undergo the pitching motion while the roll is ignored. So are 
the vibrations of the tether having a negligible mass. 

Finally, the dynamics of the two-body tethered satellite 
system with its c.m. at either of the stable Lagrange points 
namely or has been studied. Here, we explore the 
possibility of using the two-body tethered system at these stable 
points with a view to exploit the flexibility and versatility it 
offers . 


6 



CHAPTER - 2 


EXTRACTION OF KINETIC ENERGY FROM SPACE DEBRIS 


Ever since the space age began, the space has witnessed 
a rapidly escalating build-up of a large number of space objects 
such as active and spent payloads, rocket bodies and other 


mi see 1 laneous 

debris^^. Eventually, 

such 

uncontrol 1 ed 

pro 1 i f era 1 1 on of 

man-made objects is 

1 l ke 1 y 

to 

become 

a major 

source of hazard 

to the future space 

missions. 

Soon 

we must 


address to this problem before it leads to dangerous consequences. 
The possibility of a space law under the U . N . auspices warranting 
the space nations to clear one’s own debris cannot be ruled out. 
Here, we focus our attention on a more positive endeavour which 
aims at extracting significant fraction of kinetic energy of the 
wandering space objects to boost the lower payloads into higher 
orbits. In the process, the resulting reduced velocities of the 
'dead' bodies may drown them to earth-bound spiralling 
trajectories leading to their eventual demise through friction 
heating in the relatively higher density environment. 

A direct smooth transfer of energy from one body to 
another is not feasible in practice. But tether makes this 
transfer feasible due to its flexibility. Here, the basic concept 
of such a transfer has been proposed and a simple planar analysis 
is undertaken to illustrate how it works. 
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2. 1 Problem Description: 


Here, we examine the dynamics of a satellite -hook 
tethered system during the capture maneuver when the hook is made 
to attach itself to a wandering space object having a 
significantly higher speed but no practical utility whatsoever. 
At the end of such a capture maneuver results in stretching and 
tilting of the tether, which itself is later cut-off. At the end 
of the sequence of these operations the momentum transfer from the 
'dead’ body results in a significantly enhanced kinetic energy of 
the main satellite body. The planar analysis undertaken here 
mainly aims at estimating the percentage gain in the satellite 
kinetic energy. 

The available data assumed is as follows (Fig.l). 

The Satellite-hook tethered system : 

Mass = 

Semi -major axis = a^ 

Eccentricity = e^ 

The pitching angle, pitching 
rate at t = 0 

Length of the tether 

Distance between c.m. and 

Distance between c.m. and m^ 


} 


V> o , V'q • respectively 


e 



The orbiting *dead * body: 


Mass = m„ 

2 

Semi-major axis = a^ 

Eccentricity = e^ 

angle of the wandering space object as measured from the 

d satellite position at t =0, the timeof perigee passage 
given by ©£o’ 

We assume that the fast moving 'dead' body approaches 
the hook of the tethered system, and when sufficiently close, gets 
attached to it through a suitable capture maneuver that avoids 
direct hit. Subsequently, the fast moving hook - A dead' body 
combination tends to take the satellite along with it causing 
clockwise pitching of the tethered system and simultaneously 
transfering a part of the kinetic energy to the satellite in this 
process . 


The details of the geometry of motion and the 
co-ordinate axes used are shown in Fig. 2. The point *0* is the 
earth's centre while P represents the orbit perigee. The symbols 
<r,©> denote the polar co-ordinates of the c.m. of the tethered 
system at time *t' when the hook has just captured and combined 
with the 'dead' body. The corresponding length of the radius 
vector of the satellite and the attached body are represented by 
r^ and r^. The pitch angle is shown as y while ot and ft provide 
the angular separation between the two masses on either side of 
the c.m. 
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The unit vector u r is along the y-axis taken in the 

A 

direction of local vertical while u^ is along the x-axis in the 

transverse direction. The velocity vectors for the two bodies 

before and after the capture are denoted by the vectors u^ u^ and 

V 1 V 2 respectively while the hook and the tether masses are 

assumed to be relatively negligible. The symbols u u denote 

x 1 , y 1 

the components of the vector u , 1 = 1,2 along x and y axes 
respectively. Likewise v^ v represent the corresponding 
components of the velocity vector v , i = 1,2. 

2.2 Deter mi nation of post-capture velocities of masses: 

First, the principle of conservation of linear momenta 
is applied which leads to : 

m l v l + m 2 V 2 = m i u ] + m 2 u 2 = p (sa y 5 

Resolving the momentum components along x and y directions. 


X 

, V + 

1 xl 

X 2 

V x2 

P 1 

(2. 

. 1 > 

X 

, v T + 

1 y 1 

X 2 

v y2 

P 2 

(2 . 

.2) 

where X^ = 

/ (m^ 

+m 2 ) 

; X 2 


(2. 

. 2a> 


Next, we apply the law of conservation of angular momenta of the 
system about the earth’s centre, which leads to 


r, x m.v, + r_ x m.v, = r, x m.u, + r„ x m 0 u. = H (say) 

1 11 2 22 1 11 2 22 

and hence 
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-X,r, cos/? v , - X,r, sin/? v , - X„r. cosa v _ + X_r 0 sina v . = H 

11 xl 11 yl 2 2 x2 22 y2 


(2 . 3 > 


Finally, we use the energy balance. In view of the 
absence of the direct hit, the loss of kinetic energy during the 
capture maneuver would be relatively small. The main sources of 
these losses here are structural damping in the tether and some 
stress l ng-unstress l ng sequence in the vicinity of the points of 
attachment in the bodies. The actual values of losses would vary 
depending upon the materials chosen for the tether and the bodies. 
Here, we take it into account by introducing a factor K defined as 
foil ows : 


KE after the capture maneuver is completed 
KE before capture maneuver begins 


It is now easy to write the energy balance equation: 


2 2 2 2 
X, (v , + v . ) + X, (v _ + v _ ) = k. 

1 xl yl 2 x2 y2 1 


<2.4 ) 


where kj = k [ XjCu^ 2 + u yl 2 ) + ^Cu^ 2 + u y2 2 )] 

Solving the four simultaneous equations from (1) to (4) for the 


unknowns v^^ v y ^ v y ^ we obtain the following solution 


-b ± -| b 2 - 4ac 
v xl = fa 


y l 

= (k 2 \ 

'xl - «4>'K, 

> 

x2 

= Cpj 


y2 

= Cp 2 

-X 1 v y] ) / \ 2 ^ 


(2.5) 
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where 


a = t 1 + k 2 2 /k 3 2 ] 

b 2[ (Xj/xp/k^] [ k £ k 4 + Pjk^ 2 + P 2 k £ k ? ] 

c = (X ^ /X £ ) [k 4 2 /k ? 2 ] 4 ( Pl 2 4p 2 2 >/\ 2 4 2<X 1 /X 2 > 

P 2 <>W " k i 

k 2 = XjT 2 COSO( - XjTjCOS/? 

k^ = Xji^sina-XjrjSin/? 

k 4 = H 4 p 1 r 2 cosot - p 2 r 2 sinct 


2. 3 Computational Scheme : 

The steps followed are given below: 

(i> Computation of r, &, r , 9 : For this the following standard 

Keplerian relations were utilised. For a given time t, the 
mean anomaly M was obtained using 
M = n(t-t) 

l 

where 

n = Jk/a 5 

t. = time of perigee passage; tj = 0 for the tethered 

system and t 2 0 for the "dead’ body. 

K = Earth's gravitational constant. 
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The Kepler’s equation 


M = E-esinE 

was then used iteratively to solve for E to a sufficient accuracy. 
The values of r and © for the c.m. of the tethered system are now 
computed using 


tan <©/2 ) 



tan (E/2) 


and r = a(l - e cos E> 


These are next utilised to calculate r and © 


from 


f = e jK/Ca(l-e 2 )] sin© 

' I 2 2 

9 = \K a(l-e )/r 

Following the above sequence of computations, we thus obtain r 
r© for the 'dead’ body as well as the c.m. of the tethered 
system. To obtain the velocity components of the satellite, we 
would be required to take into account the librational motion as 
well which is therefore examined next. 

(ii> Planar Librational motion of the tethered system: 

The equation governing the planar librational motion^^ 
can be written as 

<1 + e cos©) ip" - 2e sin© <1 + ip'') sin© + 3 siny' cosv' = 0 
For small values of ip, the above equation can be simplified and to 
first order it can be written as 

ip" + 3 siny' cos ip - 2e sin© 
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which has the solution 


y = c^ cos (■/?©) + c^ sin <VT©> + e sin© 

y' = 3[ c^ cos(VT©> - Cj sin(VT©>] + e cos© 

For the small amplitude pitching motion we can utilize 
these approximate analytical results. Otherwise, we would be 
required to numerically solve the pitching equation for 
integration. It is now easy to obtain y as 

y - y* . jKT/r 2 

( 111 ) Computation of r^ r^ ©^ . This is accomplished using 
simple cosine law-. 


[r + 1 , - Zrl, cosy] 


1/2 


[r +1 + 2rl cosy] 


1/2 


© 


h 


© + sin 


”l 2 s l ny "I 

7 ^~ ] 


where r^ , ©^ refer to the position of the hook when the capture 
maneuver begins. 


(iv) Checking the conditions for capture maneuver: 


The conditions that must be satisfied at the point where 
capture maneuver is to be initiated are 
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and 


e 2 > e h 


[r 2 2 + r h 2 - 2r 2 r h cos<© 2 - © h >] 1/2 < € 


where e is a small length say of the order of a few meteres. It 
may be pointed out that the point P corresponding to t = 0 has to 
be so chosen that the 'dead' body is then at a location which 
eventually leads to simultaneously meeting of the above 
constraints. In the computer program we achieve this by Varying 
the time t 2 , the time of perigee passage for this body and 
eventually accepting the value of t 2 which satisfies the above 
constraints simultaneously. Here, in the absence of any specific 
knowledge of the capture mechanism, we take e = 10 m. 


< v> 


Computat i on 

For this we 


of velocity components u^j u^ 

use the following relations: 

* » • 


xl 

tO - 

1^ (0 +y )COSy 


• 

* • 

yi 

r + 

lj (0 + yOsiny 

x2 

T z e z 



• 

r 


y2 

2 



These values are computed for the instant of 
the capture maneuver occurs. 

(vi) Computation of velocity components v . v , 

a i i y * i 


u 


x2 


y2 


1 1 me 


x2 


at 


y2 : 


wh i ch 


These are obtained using the solutions given by Eqn. (2.5) 
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< V 1 1 ) 


where 


Computation of the orbital changes when the tether is 
cut-off from the satellite after the capture maneuver. 


Here, we use the following relations 


Ae . 


<e i e i } 


% increase in a 2 = [ <a’ - a^/aj] x 100 

2 2 
*■ u , j 


% gain in KE of satellite 


100 


Ae, 


6 2 6 2 


% decrease in a^, = 


£<a* - a 2 >/a 2 j x 100 


2 2 
v 2 - u 2 

% loss in KE of the 'dead' body = ( ^ |x 100 


a! 

l 


(K r >/(2K - v. r > 

l l l 


u 


2 2 j 

sin y + cos y 

f \ i | 

2 

r r v 

i_ 1 l 

fl 

L 1 « 

J 



• 


i = 1,2 


i ~1 7 2“ 

u + u . 
xi y i 


| 2 2 

v . + v 
1 x i y i 


r, 


tan ^ | v / v I 

I yi xi J 
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2.4 RESULTS AND DISCUSSION : 


The simulation results are presented in tables. Table. 1 
shows the effect of tether length on changes/gains in the orbital 

parameters. Although the use of tether connection is of vital 

importance for avoiding the direct hit and associated possibility 

of damage, the smaller tethers would be preferable for larger 

kinetic energy gains. The corresponding changes in shape and size 

of the orbit are considerable and also exhibit some dependence on 

the tether lengths. Besides, the parametric changes and kinetic 

energy gains are also significantly depend upon the orbital 

parameters of the space-wandering body [Table. 2] as is naturally 

to be expected. 

It may also be of interest to see how the trajectories 
of the floating debris changes as a result of the above momentum 
transfer. These results are shown in Table. 3. The changes in a 
and e seem to be quite substantial as also its kinetic energy 

losses. Depending upon the proposed plan of space mission and 

% 

knowledge of the various space-wandering objects, it seems 
feasible to design for a suitable capture maneuver leading to twin 
benefits with the satellite raised significantly higher altitudes 
leaving the 'dead ’ body behind to incinerate as a result of 
heating due to friction. 


17 



CHAPTER - 3 


DYNAMICS AND STABILITY OF TETHERED SPACE SYSTEM AT L AND L _ 

A 5 

The particular solutions for the three-body problem were 
determined by Lagrange in 1772. He showed that the orbit of every 
secondary revolving round its primary, has five special fixed 
points in their relative frame which are called Lagrange points 
at which the resultant of inertial and gravitational forces is 
zero. Three of these points denoted by L L and L lie on the 

*i i "■ / 

line joining the primary and the secondary. The remaining two 
points called and lying on either side of this line form an 
equilateral triangle with the primary and the secondary <Fig. 3). 
The points and are of special significance being stable. In 
other words if a small particle with suitable precise initial 
conditions is kept at either of these two points, it continues to 
stay in the vicinity of these equilibrium positions when subjected 
to small disturbances. For certain applications, the use of these 
points has been proposed in the past and its dynamics has been 
studied. Here, we propose the possibility of using a two-body 
tethered system at these Lagrange points with a view to exploit 
the flexibility and versatility it offers. For an effective 
application to be feasible, the points should be stable 

with reference to such a tethered system. This section examines 
the dynamics of a tethered system with its c.m. at the point 
and assesses its stability in the small under various parametric 
cond 1 t 1 ons . 
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3. 1 Problem Description : 


Here, we consider a two-body tethered system with c.m. 
lying at one of the 'Stable* Lagrange points, namely L 4 or L ? . 

The details of the geometry of motion and its 
co-ordinate axes used are shown in Fig. 4. For convenience 
special units are chosen, the details of which are as follows: 


(l) The distance between two massive particles and is 

taken as unity, i.e., P^P^ = 1. 

<n> The sum of masses at P J and P £ is taken as unity. Thus if 
the mass at P 2 is /j , then the mass at P^ is 1 -/u . We assume 
fj < 1/2 without loss of generality, 
lii) The unit of time is so chosen that the Universal 
Gravitational Constant G becomes unity. Let the tether 
connected masses at and be denoted respectively by 
and r> 2 which are being treated as infinitesimally small in 
compar ision to fj ^ and and which do not affect the motion 
of the massive bodies at Pj and P ^ . Here, we consider the 
particular case in which these bodies execute circular 
motion around the common system c.m. represented by '0*. 
Hence, we can easily show that 

0P 1 = h 2 = ^ 

0P 2 = ^1 = l ~V 

n = 1 


where 

n = angular velocity of the massive bodies about the system 
c.m. 
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Now, we define an inertial frame ? -77 with £ passing through 0 and 

coinciding with the line initially i.e., at t = 0 . The 

co-ordinate axes x,y denote a rotating frame with x always taken 

along the line PjP^. Then the x-axis would make an angle & = t 

with respect to £-axis at time s t' . At this general time, t, the 

tether length makes an angle yj with respect to the x-axis. The 

unit vectors along £, 77 axes are denoted by 1 , ] and that along x, 

y axes are denoted by u and u 

x y 

The meanings of other symbols used in the figure are as 

foil ows : 

= 0Q i i =1.2 

r = the radius vector of c.m. of the tethered bodies. 

r j = the radius vector P^Q, 1 = 1.2 

r 5 = the radius vector P.Q ; i = 1,2 ; 1 = 1,2 

1 1 j 

u = r /r ; i=l,2 

1 11 

C. = I /I ,i=l,2 

1 1 1 

a. = the angle between r and v ; 1 = 1,2 

1 11 

<f> = the angle between r and v^ 

For the massless tether assumed here, 

"lVl + "zVz = 0 
Vl = - v z i z 
v 2 = " V 1 
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3 * 2 Formulation of Orbital Equations of Motion for TSS : 


Let the co-ordinates of Q be <2f , -q , £) then the equation 
of motion of Q is given as follows: 


- 1 
V, T 


iv '' v ^ = ’A "■ + 

1 

After a considerable manipulation, we get 


P 2 r 


< r 2 > 

l 


] 


r A* 

r .. 

f _i 

_i I 

[ 

3 J 

v. r 


l 



3 2 

- E 

1 = 1 


2 


<p 


P /• r 5 r -v V 

' r 1 )} 

i i • /J 


On evaluating the pair of terms in the small paranthesis and 
neglecting the terms of the order of 1^ and smaller, after 
considerable manipulation we can express the above equation in the 
simplified form as 




2 2 2 
r- + ° /2> *t 


I (l-5cos 2 c] -i 
1 v J r 


r, , 2 

? 2 - 

r "1 

^2 V 

^3 - 5cos 

+ Zv i 

T 2 

— coso* x 
r l 

+ COSO* 1 

r 2 - 


<3 . 1 > 


where v 


\> 
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Now for a pure three body system without the tether or when the 
tether length is made zero, the simplified equation of orbital 
motion of the particle masses at Q itself becomes 


r 


o 


AVi 



^2 r 2 



(3.2) 
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The perturbation motion due to the tether length being finite 
and not zero is therefore obtained from 


6r = r - r = (3/2) v, 1 


2 r 

1 r i l 


1 -5cos 


S) 4 + ^ 

lJ r* r 2 


\ _ 2 T 2 

1 -5cos I — _ + 2 v 
- J « 

2 


r "i 
i 1 1* 


4 coscxj + — co so* ^ 


(5.3) 


Considering the special case when Q is nominally at L or L 

4 5 

itself when 


1 and 


60 + yj ; ex, 


120 + y 


the above perturbation equation takes the form 


2 , , „o 


6r - y r - 5 r^ cos^ ( 60°+ y >+/J 2 r 2 cos^ ( 60°- 


+ 2Vj cos (60°+ v* ) - n z cos ( 60°- y*>\ ....(3.4) 


where , 


H *•— 


C 3 . 4a > 


In order to examine the orbital perturbations, we now need to 
study the pitching librations of the tethered system which is 
considered next. 

3. 3 Equation of Pitching Motion : 

On taking the moment about the c.m. of tethered system, 
i.e., Q, the equation of pitching motion can be written as 


I W - y> ) 
zz 


However as the nominal orbit is assumed to be circular, 
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e = o 

and hence the librational equation simplifies to 


I y 

zz 


M 


zz 


where the external moment component M about Q can be evaluatec 
as 


M 


zz 


*1 X + ^2 ; 2 1 /< r 2 1>5 ] 


' 'z * “2 [ ' J 2 r l 2/<r l 2> 


3 . - 2 , , 2 . 31 

+ r 2 /(r 2 ) J 


After considerable algebraic manipulation, this expression is 
obtained in a more useful form when the equation of pitching 
motion reduces to 


v = J- <3/2>^ 1 l a l i ./<^ 1 lj+^ 2 1 2 > ] [ r sin ^ {<l-2^>cosv> -VT sinyj 


+ sin Zy 


] 


.. < 3 . 5 > 


For the motion under consideration with the c.m. of the tethered 
system at 1_ 4 or , the pitching equation can be written in a 

simplified form as 


V 


where 


sin Zy - A £ cos Zy 




. . . < 3 . 6 ) 


A = /? U-2/j)A 1 
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This is a hi gh 1 y non!, near second order ordinary differential 

equation and does not admit a known closed-form solution. So for 

a better appreciation of the system rtvnam.r-o 

ybiem oynamics, we examine the 

equilibrium positions and their stability in the small. 


3 * 3 * 1 Stability Analysis of the pitching m otion: 

First, we obtain the equilibrium configurations by 
equating the right hand side of Eqn.<5.6> to zero. This leads 
following equilibrium positions. 

V> = j tan" 1 [ ( 1 -2/j )] 

We consider the stability of the equilibrium configurations in 
some special cases: 

Case 1 : jj = fJ 1 = = 1/2 

Here the equilibrium configuration is given by 

y = O.n or ± n/Z 

< a > y = 0: 

Subst i tut i ng 

V> = u + 0 or u+tt 

leads to the following pitching equation 
u = Aj sin 2u 

Linearisation around the equilibrium position reduces it to-, 
u - 2AjU = 0 

This suggests, the system is unstable. Even a small pitching 
disturbance leads to an exponential pitching growth initially. 



( b ) yj = ± n / 2 : 

Lb! us consider the motion around yj - + tz/Z 
Subst 1 tut 1 ng 

V = u ± Tt/Z 

and linearising around u = 0 reduces the pitching equation to 
u = 2AjU 

Evidently, the equilibrium configurations y = ± tz/Z are stable. 
The frequency of persistent non-growing oscillations resulting 
from an initial distrubance is given by 

f = Cl/2n) 4 = Cl/2n> Jc| + v 2 l^> c.p.s. 

Case 2 : 

l-i £ = fj << 1 . The equilibrium pitching configurations 

are given by : 

yj = C ?t/ 6 , 7n/6), (-n/3 or 2n/3) 

We now check the stability of each of these equilibrium 
conf i gurat l ons . 

(a) y - n/6 or 7rz/6 : 

Substituting y = u + n/6 or u + 7rr/6 

and linearizing around u = 0 reduces the pitching equation 
to 

u = 4A^u 

which represents unstable pitching motion. 

<b) yi - - tt/ 3 or 2nr/3 

Substituting y - u - n/ 3 or u + 2rr/3 and linearising around 
u = 0, the pitching equation reduces to 
u* = - 4AjU 
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which implies stable/persistent non-growing oscillations with 
a frequency given by 


f = <1/2 rr) = <l/2n) J <3i>i 1 l^/<Vj 1 j + v^l^) c.p.s 

3. 4 Orbital perturbations of the tethered system : 

The orbital perturbations can be significantly 
influenced by pitching motion. However, in general, for practical 
applications, the pitch must be controlled to precise values, we 
now consider some particular cases : 

Case I Pj = /u 2 = v = 1/2 , V = ¥ c ~ 77/2 

The perturbation equation can then be simplified to 


<5 r 



< 3 .7) 


where 

v = <15/8>v,l + ■ • • ■ (5 ' 7a) 

0 T t 

and hence its components along the inertial frame can be written 
as 

< 5 ^ = - /2> sint 

Sr) = (V3v /2> cost 

C = 0 

Integration of these equat i ons 1 eads to : 
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< 5 ? 


y e [ ~ "/V2 sint + Cjt + c 2 ] 


6r> 


v e [ “ ~f\/Z cost + c^t + ] 


<$C 


C 5 t + C 6 


where c ^ , c^ , 


are constants of integration. 


Transforming the above into the relative frame xyz, the 
perturbations due to the finite tether length can be written as 


<5x = i> <c,t + c„)cost + <c,t + c,)sint 

e 1 2 3 4 


<5y = v - VT/2 - (c,t + c„)sint + <c T t + c.)cost 

e 12 ? 4 


6z 


c 5 t + c 6 


On using the initial disturbance-free perturbation condition 


x = y = z 
o o o 


X = y = z =0 
o o o 
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the orbital perturbations for the undistrubed case are obtained in 
the final form 


6x = v e -{V 3/2 sint - -/3/2 t cost} 

<5y = V e { " ~^ ,z ~ cost) + VT/Z t sint } 

6z = 0 

Case 2 = y z = y = 1/2 , y> = 0° 

The perturbation equation taken the form 


(3.8) 


(3.9) 


. . . . (3.10) 



with a new value of x> 

e 


which is now given by 



(3.11 ) 


The solution to Eqn.(3.11) can be expressed exactly as shown above 

provided this new value of x> is made use of in the solution. It 

e 

is easy to consider some degenerate cases when the mass of the 
secondary becomes negligible i.e. y = y^ = 0* 


3,5 RESULTS AND DISCUSSION : 

The above analysis clearly shows the existence of 
several stable equilibrium pitching configurations for the 
tethered systems. In the case of primary and secondary bodies 
having identical masses, the tether configuration with its length 
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perpendicular to the line joining these two bodies is stable. On 
the other hand in the situation when mass of the secondary is 
relatively negligible in comparison to that of the primary, the 
tethered system aligned with the line joining its c.m. with the 
primary represents the stable configuration as expected in this 
reduced case. However, the frequencies of oscillation in these 
two different cases differ by a factor of YI . 

In general, the orbital perturbation due to finite 

tether length causes the c.m. of the tether to drift away from the 

position of or L^. For the practical cases of interest 

considered where the pitching angle is precisely controlled to 

some specific value using a suitable control strategy, the 

perturbations present comprise of periodic as well as secular 

terms. The parameter u governing the levels of perturbation 

6 

magnitude is proportional to the square of the tether length. 

It is interesting to note that the perturbations are 
always present. Fortunately, these being only of the second order 
in dimensionless tether length (i.e. 1^/r), - which itself is 

extremely small - the orbital departures over a practical 
life-span of the TSS would be relatively insignificant. 
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CHAPTER - 4 


CLOSING REMARKS 

Some important features of the analysis and results 
based on them are summarized below. 

4.1 Kinetic Energy Transfer Analysis: 


It is posssible to exploit the kinetic energy of space - 
wandering-objects for a significant augmentation of velocity 
through the tether. Normally, small tether lengths are preferable 
for a more efficient energy transfer. It, however, requires a 
careful capture maneuver through a suitable hook. The tensions 
generated and consequent stressing of material of the tether and 
satellite have not been taken up in the present analysis. The 
analysis is restricted to in-plane motion. Evidently, there is 
considerable scope for making this investigations more consummate. 

4.2 Tethered Space System Stability at and L g : 


The pitching equilibrium configurations for the tethered 
space system kept at its c.m. at L 4 or with zero relative 
velocity have been identified and corresponding perturbation 
motion determined. In general, second order perturbations 


cause 



the TSS c.m. to continually drift away from the initial position - 
namely or L^. Since the value of the tether length in the 
dimensionless form is extremely small, the orbital departures even 
over a relatively large life-spans of the TSS are insignificant. 


CENTRA l uqpADy 
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FIG.1 TWO BODY TETHERED SPACE SYSTEM COMBINING WITH 
THE 'DEAD' BODY 
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FIG. 2 GEOMETRY OF MOTION OF A TWO BODY TETHERED 
SPACE SYSTEM AND THE 'DEAD' BODY 
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FIG. 4 GEOMETRY OF MOTION OF A TWO-BODY TETHERED SPACE 
SYSTEM AT A GENERAL POINT AROUND THE PRIMARY AND 
SECONDARY BODIES 
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Table 1 


Effect of changes in tether length on orbital parameters of 
the sate 1 1 1 te : 

a 2 = 7000 km , e 2 * 0.05 , a £ = 8000 km , e £ = 0.169 


S . NO . 

tether 

length 

change 
in e 2 (Ae^ ) 

l ncr ease 

in a^ (.%) 

gain in 

KE <*> 

1 

1 . OOOOOE+02 

1 . 09835E-01 

1 . 28450E+01 

1 .03227E + 01 

2 

2 . OOOOOE + 02 

9 . 9891 5E-02 

i 

1 . 15599E+01 

9 . 40053E+00 

3 

3 . OOOOOE+02 

i 

9 . 57170E-02 

1 . 10889E+01 

9.05963E+00 

4 

4. OOOOOE+02 

9 . 35371E-02 

1 . 08118E+01 

8.85754E+00 

5 

5 . OOOOOE+02 

9 . 40565E-02 

1 . 07302E + 0 1 

8 . 80280E+00 
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Table 2 


Effect of changes in the orbital parameters of the space— wander 1 ng 
objects on the satellite orbital parameters : 

tether length = 300 m , aj = 7000 km , e = 0.05 

rrij = 10^ kg , m^ = 10 7 kg 


S . NO . 

a 2 

6 2 

change in 
e^ (Ae^) 

i ncrease 

in a^ (%) 

gain in 

KE <%) 

1 

8 . OOOOOE+06 

1 . 69000E-01 

9.89754E-02 

1 . 1 3278E+0 1 

9 - 24490E + 00 

2 

8 . 10000E + 06 

1 . 79000E-01 

1 . 05604E-01 

1 . 18142E+01 

9.6201 5E+00 

3 

8 . 20000E+06 

1 . 89000E-01 

9 . 39943E-02 

1 . 0 570 6E+0 1 

8 . 68642E+00 

4 

8 . 30000E+06 

1 . 99000E-01 

1 .071 3 5E-01 

1 . 18509E+01 

9 . 64536E+00 

! 

5 

8 . 40000E+06 

2 . 09000E-01 

1 . 09883E-01 

1 .20178E+01 

9 - 77202E+00 

6 

8 . 50000E+06 

2 . 1 9000E-01 

1 . 08786E-01 

1 . 14976E+01 

9-401 91E+00 

7 

8 . 60000E+06 

2 . 29000E-01 

1 . 16054E-01 

1 .236B0E+01 

1 . 00342E+0 1 

8 

8 . 70000E+0 6 

2 . 39000E-01 

1 . 18073E-01 

1 . 24635E+01 

1 . 01 053E+01 

9 

8 . 80000E+06 

2 . 49G00E-0 1 

1 . 1 5547E-01 

1 . 18100E + 01 

9 . 63493E + 00 

10 

8 . 90000E+06 

2 . 59000E-01 

1 .22217E-01 

1 . 26556E + 01 

1 . 02458E+01 
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Table 3 


Effect of changes in tether length on orbital parameters of the 
space— wander i ng objects : 

a 2 = 7000 km , e ][ = 0.05 , a £ = 8000 km , e £ = 0.169 


S . NO . 

tether 

length 

change 

in e 2 

decrease 

in <^> 

loss in 

KE <30 

1 

1 .00000E + 02 

5.82293E-04 

8.021 87E-02 

9.5751 5E-02 

2 

2 . 00000E+02 

2 . 32965E-04 

8 . 90625E-02 

8.73239E-02 

3 

3. 00000E + 02 

1 . 451 37E-05 

8 . 46875E-02 

8 . 41 92 4E -02 

4 

4. 00000E + 02 

-5 - 91 427E-04 

1 .06850E-01 

8 . 22626E-02 

5 

5 . 00000E+02 

-6 . 60300E-04 

1 . 07056E-01 

8 . 1 6756E-0 2 
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